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Abstract 

The utility-based pricing of defaultable bonds in the case of stochas- 
tic intensity models of default risk is discussed. The Hamilton- Jacobi- 
Bellman (HJB) equations for the value functions is derived. A finite 
difference method is used to solve this problem. The yield-spreads 
for both buyer and seller are extracted. The behaviour of the spread 
curve given the default intensity is analyzed. Finally the impacts of 
the risk aversion and the correlation coefficient are discussed. 
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1 Introduction 



Credit risk models can be partitioned into two groups known as structural 
models and reduced form models. 

Structural models were pioneered by lMertonl (11974I ). The basic idea, common 
to all structural-type models, is that a company defaults on its debt if the 
value of the assets of the company falls below a certain default point. Such 
a default can be expected. In these models it has been demonstrated that 
default can be modelled as an option and, as a result, researchers were able 
to apply the same principles used for option pricing to the valuation of risky 
corporate securities. The application of option pricing avoids the use of risk 
premium and tries to use other marketable securities to price the option. The 
use of option pricing theory set forth by Merton offers much more accurate 
prices, but provides information about how to hedge out the default risk. The 
shortcoming of the structural approach is that it underestimates credit risk 
in the sense that the corporate bond is cheaper than the default free bond 
ev en if the f irm v alue is larger than the default point. Subsequ ent to the work 
of iMertonl (119741), there have been m any extensions such as iBlack and Cox 
(119761 ) . iLongstaff and Schwartzl (119951 ). etc. The second group of credit mod- 
els, known as reduced form mod el s, are more recent. Th e se mo d els, including 
the lArtzner and Delbaenl (119951) . Jarrow and Turnbulll (119951 ). iJarrow et al 



( 19971 ) and iDuffie and Sing letonl (I1999h models, do not look inside the firm 
Instead, they assume that default occurs without warning at an exogenous 
default rate, or intensity. So, the default event would be unexpected. The 
dynamics of the intensity are specified under the pricing probability. In- 
stead of asking why the firm defaults, the intensity model is calibrated using 
the market data such that the credit rating, the stock prices, risk free rate 
and etc. As the market data come not only from credit risk but also from 
other factors, thus the reduced form model overrates the credit risk. This 
is the shortcoming of the reduced approach. One of the recent trends i s 
to combine bot h ap proaches, see for examp le ICathcart and El-Jahell (120031 ). 
Realdonl (|2006[ ). In ICathcart and El-Jahell (|2003|) a PDE method is used, 



providing a semi-analytical pricing formula of defaultable bond combining 
the two approaches when the short rate follows CIR model and the default 
intensity is linearly dependent on the short rate. The case with exogenous 
default recovery and no correlation between the firm value and the short rate 
is studied. It is assu med that the ex pected and unexpected default recovery 
rates are the same. iRealdonl ( 120061 ). using PDE method, provided an ana- 
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lytical pricing formula of defaultable bond with both approaches in the case 
that the default intensity and the short rate are constant and the default 
recovery is endogenous. He also provided an analytical pricing formula when 
the short rate is uncorrelated with the firm value and the default recovery 
is exogenous. He didn't assumed that the unexpected and expected default 
recovery rates are the same. 

However, a market with stochastic default intensity is incomplete in that the 
intensity is a source of uncertainty that is not traded. The incompleteness of 
the market usually gives rise to infinitely many martingales measures, each 
of which produces a no-arbitrage pricing. So, it is not clear which one to 
use in the pricing of the defaultable bond. The superreplication price (upper 
hedging price) is the minimal initial wealth needed for hedging the claim 
without risk. If the bank decides to charge a superreplication price for selling 
a bond so that it can trade to eliminate all risks, the price is usually high and 
unrealistic. There are three major approaches that have been developped in 
searching for solutions of pricing and hedging in incomplete markets. 
One is to pick a specific martingale measure for pricing according to some 
optimal criterion. For example o n can cite the minimal relative entro py mar- 
tingal e measure in lAvellanedal ( 119981 1 or the q-optimal measure in iHobson 
(120041 1 . The drawback is this approach sometimes provides a result that is 
not very reasonable financially. Another approach is to identify the martin- 
gale measure consistent with the market price. The problem is it does not 
provide a hedging strategy at all. Instead of choosing an equivalent mar- 
tingale measure for valuation, a dynamic utility-based valuation theory has 

been developed producing the so-called Indifference price. 

Utility indifference pricing was first introduced by iHodges and Neuberger 
(119891 ) . It is an alternative where the price is uniquely determined at the 
cost of depending on the preferences of the pricer. The writer utility indif- 
ference price is the value of the initial payment that makes the seller indif- 
ferent to whether to sell the contract or not. From practical point of view 
utility indifference pricing has at least two advantages: 1) it does not refer 
to the market portfolio; 2) it generates an optimal hedging strategy in the 
sense that the resulting utility is equal to that of an optimal pure invest- 
ment. The mathematical structure of utility in difference pricing h a s bee n 



Delbaen et al 



well character i zed by numerous resear c hers, cf. iRouge and Karoui ( 2000 ). 

" (120021) iBechereri (l2003h . iMusiela and Zariphopouioul (l2004h . 



Mania and Schweizerl (120051). and references therein. 



Bielecki et al.l ( 120041 ). ICollin-Dufresne and Hugonnierl ( 120011 ). ISircar and Zariphopoulou 
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(120071 ) , IShoudal (120061) and ILeung et al.l ( 120081 ) apply the utility-based pric- 
ing to credit risk. iBielecki et al.l (120041 ) discussed the utility-indifference price 
of defaultable claims within the reduced approach by a backward stochastic 
differential equation. They studied a particular indifference price, based on 
the quadratic criterio n and solved the problem by th e dual ity approach for 
exponential utilities. CoUin-Dufresne and Hugonnierl ( 200ll ) study the util- 
ity based prices of Event Sensitive Contingent (ESCC) claims under two 
scenarios of resolution of uncertainty for event risk: when the event is con- 
tinuously monitored or when it is revealed only at the payment date. In 
both cases, they transformed the incomplete market optimal portfolio choice 
problem of an agent endowed with an ESCC into a complete market problem 
with a state and possibly path dependent utility function. They also obtain 
an explicit representation of the utility based prices under both informa- 
tion resolution scenarios in the c ase of negative exponential utility function. 
Sircar and Zariphopouloul ( 120071 ) apply the utility indifference valuation in 
intensity-based models of default risk where the default time is the first jump 
of a time-changed Poisson process. They derive the Hamilton- Jacobi-Bellman 
(HJB) equations and analyze resulting yield spreads for single-name default- 
able bonds, and a simple rep resentative two -name credit derivative when the 
default intensity is constant. IShoudal ( 120061 ) discussed the utility-based pric- 
ing of defaultable bonds where their recovery values are unpredictable. He 
derived a partial integro-differential equations that the utility-based bond 
price solves. He also extracted credit risk premium from the yield spread of 
defaultable b onds and classified them to default-timing risk, recovery risk and 
spread risk. ILeung et al.l ( 120081 ) apply the technology of utility-indifference 



valuation for defaultable bonds in a structural model for Black-Cox-type. 
They derive the HJB equations and simplify them to the linear (Feynman- 
Kac) differential equations. Finally, they find that the utility valuation has 
a significant impact on the bond prices and yield spreads. 
The indifference valuation should be attractive to participants working in the 
OTC market. It is a direct way for them to quantify the default risks they 
face in a portfolio of complex instruments, when calibration data is scarce- 
Motiv ated by theses issues, ou r paper is similar to ISircar and Zariphopoulou 
( 120071 ) and ILeung et al.l (120081 ) but here the default time r is the first jump 
time of a Cox process with a random intensity A > 0, which is correlated 
with the firm's stock price S. However, we consider the full filtration F = 
(J r t)t>o to make the default time r a F-stopping time. We also keep the 
same notation as ISircar and Zariphopouloul (120071 ) and ILeung et al.l (120081 ) 



to facilitate the comparison of the results. Assuming the intensity follows a 
Cox et al. ( 19851 ) (CIR) model, a finite difference method is used to solve the 
nonlinear reaction-diffusion equations satisfied by the value functions. The 
analysis of the bid and ask spread curves clearly shows the nonzero short- 
term yield-spread reflectin g that the de f ault r isk comes as a surprise. This 
effect is not significative in iLeung et al.l (120081 ) even if their results enhance 
short term yield spreads compared to the standard Black-Cox valuation. 
An interesting question is now, how the yield-spread behaves depending on 
different factors. It is quite obvious, that when the default intensity increases, 
the spread curves will be high. But what kind of behavior for the spreads 
would be expected with varying maturity to the exposure, given a default 
intensity? Some empirical studies have focused on a split behavior o f the 
credit spread s accor di ng to the cr edit quality of the firm; see iFonsl (119941 ) and 



Jarrow et al.l (119971 ). IFonsl (119941 ) using the rating information produced by 



Moody's, discovers an almost strictly upward slope for the spread curves for 
bonds in investment-grade classes. The so-called humped-shaped behavior 
can be observed for th e rating BB and a do wnward slope has been calculated 
for the rating class B. iJarrow et al.l ( 119971 ) found for the issues rated AAA, 
AA or A a strictly upward sloping credit spread curve. For the issues rated 
BBB, BB and B a humped shape curve is obtained, while for the rating class 
CCC a downward sloping curve is observed. It is important to notify that in 
the precedent studies of split behavior, markets are assumed to be complete 
and free of arbitrage opportunities. 

We apply the model to analyze the shape of the spread curve in the incom- 
plete market according to the level of the default intensity. We found for the 
bid spread curve, a upward sloping curve for < A < 0.06, a S shape curve 
for 0.06 < A < 0.15 and a downward sloping curve for 0.15 < A < 1. For 
the ask spread curves, a upward sloping curve is observed for < A < 0.25 
and a humped sloping curve for 0.25 < A < 1. We also observed the positive 
relationship between the bid spread and the risk aversion parameter and a 
negative relationship between the ask spread and the risk aversion whatever 
the maturity. The impact of the correlation coefficient on the spread curves 
is analyzed. It shows the correlation coefficient increases the bid and ask 
spread for longer maturity. Finally we recover the classical structural models 
by proving they are the limit cases of our model. 

The remainder of this paper is organized as follows. Section [2] provides a 
brief introduction to utility maximization and indifference pricing. Section 
E] examines the indifference pricing of defaultable. Using the exponential 
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utility function, the Hamilton-Jacobi-Bellman (HJB) equations are presented 
and simplified to some reaction-diffusion equations. The dual optimization 
problem is revised and then some numerical results are presented. 
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2 Utility Maximization and Indifference Pric- 
ing 

2.1 Utility Maximization 

Utility maximization or the maximization of expected utility is a very basic 
topic in economics. A rigorous treatment of this subject within the fra mework 



of ma them atical finance was initiated by the two seminal papers by iMerton 



fll969f ) and ISamuelsor] fll969h . Since then the problem of utility maximiza- 



tion has found its way into numerous textbooks : from the more economic 



point of view, e.g.. lHuang and Litzenbengerl (119881) . a s we ll as from the more 



mathe matical side, e.g., iKaratzas and Shrevd (|1998) or iFollmer and Shied 



fl2002h 



Let us start by giving an introduction to the kind of questions we are going 
to treat in the following and describing the setting in which this will be done. 
We assume that there is an economic market described by a probability space 
(Q, {TtjoKtKT-, I )- & denotes the possible states of the world, J-" t the filtration 
which contains all information up time t, and P is the probability measure 
describing the likelihood of certain events. In the market there are a certain 
number of assets (Sl)o< t <T- We take one asset as numeraire, meaning that we 
assume that this special asset, called the bond or cash account, is constant 
equal to 1 and that the other assets are denoted in units of the bond. The 
other assets, in the following called stocks, are assumed to be risky and 
are described by an adapted, i.e., JF t - measurable, stochastic processes. The 
market we consider is furthermore assumed to be frictionless. This means 
that the agents who trade in the market, i.e., buy and sell assets, do not 
face any transaction costs. Thus the gain or loss of an agent engaged in 
the market following a trading strategy 9 is given by the stochastic integral 
J 9dS, also meaning that the trading strategy is self-financing, there is no 
money exogenously added or withdrawn. The agents face no constraints like 
short-selling constraints prohibiting them from holding a negative amount 
of a certain stock. In addition we always assume that the market is free of 
arbitrage, basically meaning that it is not possible to make money out of 
nothing. 

We now consider a special agent who has sold at time a contingent claim 
C that matures at time T. Such contingent claim is modelled as an Tt- 
measurable random variable. The agent will have to pay back the random 
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amount of C units of the bond at final time T. Let us also assume that his 
current wealth level is c. Thus the agent, if he chooses the trading strategy 
9, at final time T will hold 



Z{B){u) = c+ 6(u)dS(u) - C(u) 
Jo 

We call the random variable Z(9) the outcome of the trading strategy 9. 
Now, the question is how that the agent choose the trading strategy? We 
assume that his preferences h ave an expected utility representation, follow- 
ing the theory introduced by Ivon Neumann and Morgensternl (119531 ). This 
means that he tries to maximize his expected utility over all admissible trad- 
ing strategies 9: 



E(U(Z(9))) -»■ maxfl 



The Ivon Neumann and Morgensternl (119531 ) utility function U(x) is assumed 
to be strictly increasing corresponding that investors prefer more wealth to 
less and strictly concave. The concavity reflects the issue of risk aversion: 
people are assumed to dislike risk and thus prefer for example to get 100 
currency unit for sure compared to having a 50-50 chance of getting 200 
currency unit or nothing. The concavity implies : 



\U{x) + (1 - X)U{y) < U{Xx + (1 - X)y), VA e (0, 1). 
In our example |Z7(200) + |C/(0) < 17(100). 

Let us introduce the coefficient of absolute risk aversion defined by 



Ra(x) 



U"{x) 
'U'(x) 



Since U (x) is assumed to be strictly increasing, which means (provided U(x) 
is sufficiently differentiable, which we assume) U'(x) > 0, and strictly con- 
cave, i.e., U"(x) < 0, the risk aversion R Jx) is greate r tha n 0. T h e mea sure 
R a {x) of risk aversion was introduced by lArrowl (Il97ll ) and iPrattl (11964I ). As 
indicated above it is very reasonable to use the curvature i.e. the degree of 
concavity U"(x) as a measure for the risk aversion. The reason why U"(x) 
is divided by U'(x), which correspond s to a normalization, is the fact that 
von Neumann and Morgensternl (119531 ) utility functions are only unique up 
to positive linear transformations. That is U(x) and U = aU(x) + b for 
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a, h e M, a > correspond to the same preference relation of the agent. Di- 
viding U"(x) by the derivative U'(x) respectively JJ'ix) ensures the same risk 
aversion coefficient. There exist many utility functions and the most useful in 
utility maximization are power utility functions U(x) = fz^-; R ^ I, R > 
and exponential utility U(x) = — e~ ix ; 7 > 0. The description of both of 
these functions is given by describing the HARA and Non-HARA utilities. 

A utility function is of the HARA class if R a (x) = A + Bx x E Id where Id is 
the interval on which U is defined and B is a non-negative constant. 

• If B > and fi^l, this leads under some assumptions to the standard 
power utility function U(x) = The power utility has constant 
relative risk aversion (CRRA) of R, where the relative risk aversion 
R r (x) = xR a {x). The power utility has fundamental assumption that 
it requires wealth to be non-negative. 

• If B — 1, we have under some assumptions the logarithmic utility 
function U(x) = \n(x) which also requires a non-negative wealth. 

• If B — 0, this leads to the exponential utility function U(x) = —e~ lx 
which has a constant absolute risk aversion (CARA) of 7. Bigger 7 
corresponds to higher degree of risk aversion. In particular, 7—7-00 
indicates absolute risk aversion while 7—7-0 corresponds to risk neu- 
trality. In contrast to the power utility function, the exponential utility 
function allows for negative wealth. 

The Non-HARA class regroups all utility functions which do not fit the 
HARA class. An utility function which belongs to the Non-HARA is the 
quadratic utility function U(x) = x — bx 2 where b > 0. This utility func- 
tion decreases over part of the range, violating the assumption that in- 
vestors desire more wealth, but has excellent tractability properties; see 



Henderson and Hobsonl (120041 1 . 



2.2 Utility Indifference Prices 

The utility indifference buy (or bid) price p b is the price at which the investor 
is indifferent, in terms of maximum expected utility, between paying nothing 
and not having the claim Ct and paying p b now to receive k > units of 
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the claim Ct at time T. Assume that the investor initially has wealth x and 
zero endowment of the claim. Mathematically, this expressed as 

sup XTeA{x) E[U(X T )} = sup XTeA(x _ pb{k)) E[U(X T + kC T )}, (1) 

where A(x) is the set of all wealth X? which can be generated from initial for- 
tune x. Let M(x) = supx T£A ^E[U(X T )] and H b (x, k) = sup XTeA f x \E[U(X T + 
kCr)]] the utility indifference buy price p b is the solution to 

M{x) = H b (x-p b {k),k), (2) 

that is the investor is willing to pay at most the amount p b (k) today for k 
units of the claim Ct at time T. Similarly the utility indifference sell (or 
ask) price p s (k) is the smallest amount the investor is willing to accept in 
order to sell k units of Ct- The utility indifference sell price p s (k) solves the 
problem 

M(x) = H s (x + p s (k) } k) } (3) 
where H s (x, k) = sup XTeA{x) E[U(X T - kC T )}- 

From equations <^ and to compute the utility indifference price of a 
claim, two stochastic control problems must be solved. The first one is the 
optimal investment problem w hen the i nvest or has a zero position in the 
claim; this is equivalent to the Merton ( 197ll ) optimal investment problem 



and the value function is denoted by M(x). The second one is the optimal 
investment problem when the investor has bought or sold k units of the claim 
and is denoted by H(x, k). 



Utility indifference prices have some appealing properties. Theses properties 
are given in the following 

Theorem 1. 1. Non-linear pricing. Firstly, in contrast to the arbitrage- 
free prices, utility indifference prices are non-linear in the number of 
claims, k. This means that the buyer is not willing to pay twice as 
much for twice as many claims, but requires a reduction in this price to 
take an additional risk. Alternatively, a seller requires more than twice 
the price for taking on twice the risk. This property can be seen by the 
value function H since the utitlity function U is concave. 

2. Recovery of complete market price. If the market is complete or the 
claim is perfectly replicable, the utility indifference price is equivalent 
to the complete market price for k units. 
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3. Monotonicity. Let p 1 and p 2 be the utility indifference prices for one 
unit of payoff C\ and C T respectively and C\ < C T , then p 1 < p 2 . 

4- Concavity. Let p\ be the utility indifference bid price for the claim 
XC T + (1 - X)C T where A G [0, 1], then 

Px > Xp 1 + (1 - X)p 2 



We br iefly recall the proof of theses results. (See also iHenderson and Hobson 



( 120041 )) 



Proof. ([2]) Let Ct the contingent claim which is replicable, p° the arbitrage- 
free price of the claim and p b (k) the buyer indifference price of k units of 
the contingent claim. We denote by Rt the time T value of unit of currency 
invested at time 0. Since Xt G A(x), Xt = xRt + Xt for some Xt G A(0), 
where A(0) is the set of claims which can be replicated by zero initial wealth. 

*** C 

Since Ct is replicable from an initial fortune p c , we have Ct = p c Rt + Xt 
where Xt G A(0). Then for X T G A(x), X T + kC T = (x + kp c )R T + 
X T + kX T C = (x + kp c )R T + X T ' where X T ' = X T + kX T ° e A(0). Thus 
X T + kC T G A(x + kp c ). But H b (x, k) = sup XTeA(x) E[U(X T + kC T )\ = 
sw Px T eA(x+kp c )E[U(XT)] = M(x + kp c ) and thus p b (k) = kp c . It means that 
the buyer indifference price is k units times the arbitrage-free price. □ 

Proof. © Let assume that C\ < C T , then H b (x, 1, C8j) < H b (x, 1, C|), with 
H b (x, 1, C T ) and H b (x, 1, C T ) are buyer's value functions for one unit of the 



on 



claims C T and C\ respectively. As H b (x — p 1 , 1,C T ) = H b (x — p 2 , 1, C 
deduce M(x + p 1 ) < M(x + p 2 ) and thus p 1 < p 2 . □ 

Proof. (Tj0) Let X T be the optimal target wealth for an individual with initial 
wealth x — p % due to receive the claim C T - Then, M(x) = H b (x — p\ 1, C T ) = 
snp XTeA{x _ pl) E[U(X T + C T )\ = E[U{X T + (%)]. Define X T = XX T + (1 - 
X)X T , we have X T G A = A(x - Xp 1 - (1 - X)p 2 ). So, 
H b (x - Xp 1 - (1 - X)p 2 , 1, XC T + (1 - A)C|) 

= sup XT6 ^[C/(X r + ACi + (l-A)C|)] 

> + ACi + (l- X)C T )) 

= E[U(X(X T + C T ) + (1 - X)(X 2 + C 2 )] 

> XE[U(X T + C T )\ + (1 - X)E[U(X 2 + C T )} 
= M(x) 

= H b (x-p x ,l,XC T + (l- X)C 2 ) 
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Thus H\x-Xp l -(l-A)p 2 ,l,ACf + (l-A)C|) > H\x -p x , 1, AC* + (1 - 
A)C|) and therefore p A > Ap 1 + (1 — A)p 2 □ 

If we consider ask price rather bid price then p\ is convex rather than concave. 
These value functions are the solutions of the Hamilton- Jacobi-Bellman (HJB) 
equations. After solving theses equations, we can find the utility indifference 
bid price and ask price for the claim. In the next chapter, we will derive the 
utility indifference for defaultable bonds in the case of stochastic intensity. 



3 Indifference Pricing For Defaultable Bonds 

3.1 Maximal Expected Utility Problem 

Let us start with a defaultable bond of a firm with expiration date T < oo. 
Unlike in a traditional structural approach, default time r is the first jump 
of a Cox process with a random intensity A > 0, which is correlated with the 
firm's stock price S. The price of the asset S is described by a geometric 
Brownian motion. The intensity process is X(Y t ) where A(.) is a non-negative, 
locally Lipschiptz, smooth and bounded function. The dynamics of S and Y 
are 

f dS t = fiS t dt + aS t dW t \ S = S>0, , . 

\ dY t = b(Y t )dt + a(Y t ) {pdW} + p'dW 2 ) , Y = X, U 

with p' = a/ 1 — p 2 . The processes W l and W 2 are two independent Brownian 
motions defined on a probability space (Q, J 7 , P) and Tt is the large filtration 
enough to support a counting process N t = l{ T <t} and stochastic processes 
((S u , Y u ); < u < t) which are independent of N. Hence we may write Tt = 
Qt V U u where G t = a {(S u , Y u ) : < u < t} and H t = a {N u : < u < t}. 
The background filtration Q t contains all economic informations except for 
the default times available to investors up to time t. However, the default 
time r is not a stopping time with respect to Q t but it is under Tt- In addition, 
if we condition on the filtration Q t , N t is an inhomogeneous Poisson process 
with the time-varying intensity A. The parameter p G (—1, 1) measures the 
instantaneous correlation between shocks to the stock price S and shocks 
to the intensity-driving process Y. We also assume that there exists a unit 
exponential random variable E\, independent of the Brownian motions. The 
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default time r of the firm is defined by 

t = inf js > : jf X u du > #1 j . 
We make the following assumption: 

Assumption 1. TTie coefficients a(.) and o(.) are suc/i i/iat (0) /ias a unique 
strong solution Y t which lies in M. + for all t G [0, T], P-a.s. 

Observe that Assumption 1 is weaker than the usual global Lipschitz conti- 
nuity and growth conditions that are sufficient for the existence of a strong 
solution Y t . This avoids ruling out from the stochastic intensity model such 
that the CIR model, for which Y t is square root diffusion that does not satisfy 
the global Lipschipz condition. 

Assume that a bank account is also traded in the market and the risk free 
rate interest rate is r. The investor's control process is nt, the amount held in 
the stock at time t, until tAT. The control process ir is admissible i.e. it is Tt 
measurable and satisfies the integrability contraint E{j Q iigds} < oo. The 
set of all admissible strategies is denoted by A. In t < t AT, the investor's 
wealth process X t follows 

dX t = n t ^ + r(X t -it t )dt (5) 
= (rXt + TTtifi-r^dt + andWt 1 . (6) 

If the default event occurs before T, the investor has to liquidate holdings 
on the stock and deposit in the bank account. For simplicity, we assume he 
receives full pre-default market value on his stock holdings on liquidation. 
Therefore, given that r < T, for r < t < T, we have X t = X T e r( - t ~ T \ 
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Figure 1: Evolution of Wealth before and after the default event 



Using the exponential utility of discounted wealth, we are first interested in 
the optimal investment problem up to time T of the investor who does not 
hold the defaultable bond. At time zero, the maximum expected utility of 
discount payoff takes the form: 

sup^E (-e~^- rTx -h {T>T} + (-e-^-^))l {r < r} ) (7) 

To simplify the discount variable, let us apply the Ito Lemma to e~ rt X t 

d(e- rt X t ) = e- rt (7r t (ii - r)dt + <rir t dW}), 

and switching the discount variable X t — > e~ rt X t , the excess return fj, — > /j,— r 
and iT t — > Tv t e~ rt , we get 

dX t = n t (fAdt + adWl) . (8) 

Next, we define the stochastic control problem initiated at time t < T, and 
define the default time 

r f = inf |s > t : J \ u du > E^ . 

The Merton value function becomes 

M(t,x,y) = sup„ teA E t (-e-^l {Tt>T} + (- e ^)l {rt < T} ) , (9) 
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where (x,y) = (X t ,Y t ), E t (-) = E{-\T t ). 

We are now interested to the same problem of the investor who owns the 
defaultable bond of the firm, which pays 1 currency unit on date T if the 
firm has survived till then. The bond's holder value function is: 

H%x,y) = su VnteA E t (- e -^ x ^h {Tt>T} + (-e^)l { ^ n ) , (10) 
where c = e~ rT . In the same way, the bond's seller value function is : 

H'(t,x,y) = su P7VteA E t (- e -^ x -- c h {Tt>T} + (-e^)l {Tt < T} ) . (11) 
Remark 1. 

The value functio ns M, H b , H s are concave , incre asing in x and uniformly 
bounded in y . See lSircar and Zariphopoulou 



In the next section, we derive the Hamilton- Jacobi-Bellman equations satis- 
fied by these value functions. 

3.2 The Hamilton-Jacobi-Bellman Equations 

Proposition 1. Consider the Merton value function 

M(t,x,y) = su Pnt£A E t (_ e -7*r 1{7i>T} + (-e-^)l {rt < T }) , 
subject to the wealth and state variable contraints 

dX t = 7v t (fidt + adWl) 

and 

dY t = b(Y t )dt + a(Y t ) [pdW] + p'dW?) 

with the coefficients a(.) and b(.) verifying the assumption 1. 
The Merton value function M is the unique viscosity solution in the class of 
functions that are concave and increasing in x, and uniformly bounded in y 
of the Hamilton-Jacobi-Bellman (HJB) equation 

( E,m^v)n =0 (12) 
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Proof. Assuming that an optimal amount exists, we derive from 



sup nt E t (- 



-yX 7 



l{r t >T } + (-e-^)l{n<T } ) - M{t,x,y) = 



for some small h > 0, 

sup^ [E t+h (-e^ x n {rt>r} + (-e-^)l{n<T})] ~ M (M,Z/) = °- 

The expression i£t+h (— e~ 7XT l{ Tt> r} + (— e~ 7Xrt )l{r t <T}) is nothing else than 
the expected utility for an investor starting with wealth X t +h at time t + h. 
Therefore, for any control ir s with s > t + h, 

E t+h (-e^ XT l {Tt>T} + (- e ^)l{r t <n) < M (* + h,X t+h ,Y t+h ), 

with equality holds for the optimal n*. Hence, 

sup*„t<.<t+hEt [M (t + h, X t+h , Y t+h )} - M(t, x, y) = 0. 

Assuming an optimal behavior from t + h on T, the optimal n* has only to 
be determined until t + h, and not on the whole time horizon. Dividing by 
h and applying the limit h — > 0, the above equation becomes 



sup^ t Uim h ^ E t - [M {t + h, Xt+n, Y t+ h) - M(t, x, y)]j = 
Using the theorem of bounded convergence, this equation becomes 



sup, t j^lim^o- [M(t + h,X t +h,Y t + h ) - M(t,x,y)]j = 0, 
but lim^ox [M (t + h, X t+h , Y t+h ) - M(t, x, y)} = |M(t, x, y) and the value 



function M verifies 



E t [dM(t,x,y)] 
dt 



sec 



Sennewaldl (120071 ) or lMalliaris and Brockl (I1982T) for more details). 



Follow ing the arguments of Theorem 4.1 and 4.2 in lDuffie and Zariphopoulou 
( 119931 ). we deduce that the value function M is the unique viscosity solution 
of the HJB □ 
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Let us now derive the explicit form of the HJB equations satisfied by the 
Merton function M. We have 



E t (dM(t,x,y)) = 

\M t + tt/iM x . + b(y)M y + ^a 2 n 2 M xx + ^a(y) 2 M yy + pana(y)M xy ^j dt+ 

X(y) (-e-T* - M) dt, 
Dividing by dt and applying the supremum, we obtain 

M t + C y M + sup, (^a 2 7i 2 M xx + vr (paa(y)M xy + pM x )^j 

+X(y) (-e~^ - M) = (13) 
with M(T, x, y) = -e-T* and C y = \a{y) 2 j^ 2 + b{y)f y 

Similarly, the value function H b of the investor, who owns the defaultable 
bond, satisfies the HJB equation is: 

H b + C y H b + sup, (VV#1 + tt (paa(y)H b xy + pH b )^j 

+A(y) (-e-^ - H b ) = (14) 

with H b (T,x,y) = -e~^ x+c \ 

The value function H s of the seller of the bond also satisfies : 

H' t + C y H s + sup, Q^vr //;.,. + 7r {paai^H^ + /iF*)) 

+A(y) (-e-^ - = 0, (15) 

with H s (T,x,y) = -e~^ x ~ c \ 

Using (jSJ) and (j3J), the buyer's indifference price Pq(T) and the seller's indif- 
ference price Pq{T) (at time 0) of the defaultable bond with maturity T are 
defined by 

M(0,x,y) = H b (0,x-p b ,y) (16) 

and 

M(0,x,y) = H s (0,x + p s Q ,y). (17) 
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Remark 2. 



It is well-known that the indifference price under exponential utility does not 
depend on the investor's initial wealth x, which is an attractive feature of 
using this utility function. 

In the next section we factorize the value functions M, H b , H s , that leads 
to some reaction-diffusion equations. 



3.3 Reaction-Diffusion Equations 

In order to simplify the above HJB equations let us introduce the familiar 
distortion scaling: 

M(t,x,y) = -e^ x u{t,y)^?, 
H b (fx,y) = -e-^w b (t,y)^, 

and 

H s (T,x,y) = -e-^w s (t,y)^, 
with u,w b ,w s : [0,T] x R — > R + . For the proof of the factorization, see 



Sircar and Zariphopouloul (120051 ). Plugging these expressions in the HJB 



equations, we get the following reaction-diffusion equations: 

u t + Kyu - (1 - p 2 ) + A(y)) u + (1 - p 2 ) \{y)u~ e = 0, (18) 
with the final condition u (T, y) = 1. 

w b t + K y w b - (1 - p 2 ) (j^ + KV)) w b +(l- p 2 ) e' c X(y)(w b )- s = 0, (19) 
with the final condition w b (T, y) = 1. 

w} + K y w s - (1 - p 2 ) (JL + \(y)j w s +(l- p 2 ) e-^X(y)(wT e = (20) 
with the final condition w s {T,y) = 1, where 

1 — p 2 
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and 

PH . . d 
a dy 

With the simplified form of the value functions, the buyer indifference's price 
is given in the following way; from section 12. 2\ 



M(0,x,y) = H b (0,x-p b ,y) 
-^(0,2/)^ = e- 7(x - p o +c V(0,y)^ 



this yields to 



7 (1 -p 2 ) & W0, 2/) 
It is shown in ISircar and Zariphopoulou fl2007t ) that w s (t,y) < u(t,y) < 



w b (t, y) for (t, y) G [0, T] x E, so p b (T) < e~ rT . The yield spread of the bond 
is defined by 

V b o(T) = -iog(^(T)) - r, (22) 
which is non-negative for all T. 

In the same way from section 13. 2\ the seller indifference's price is 

Pl(T) = - .^log ( ^M) • (23) 



7(1 ~ P 2 ) \w s (0,y) 
and the seller's yield spread is non-negative for all T. 

The indifference price in this model can be characterized as the solution of 
the non-linear reaction-diffusion equations related to the functions u, w b and 
w s . We present in the next section the corresponding dual problems which 
is an alternative approach to evaluate the indifference price. 

3.4 Dual Problem: Relative entropy Minimization 

We recall the primal optimal investment problem for the buyer at time zero, 

V = su VneA E {-e^ x - +c h {T>T} + (-e~^)l {T < T} ) . 
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Under the framework that the underlying processes are locally bounde d semi- 
martingales and the utility function is exponential, iDelbaen et al.l (120021 ) 
proved the duality relation 

where x is the initial wealth and T> is the value of the dual optimization 
problem 

V = inf QeP/ (V {c} + Ui{Q\P)\ (24) 
where 1-L(Q\P) is the relative entropy between Q and P, namely, 



H f Q \P ) = I (§)}>Q« P > (25) 

| oo, otherwise 

where Q « P denotes Q is absolutely continuous with respect P. In fl24|) . 
Ff denotes the set of absolutely continuous local martingale measures with 
finite relative entropy with respect to P. 

Assumption 2. There is an equivalent local martingale measure with finite 
relative entropy. Thus, 

P/(P)nP c ^ (26) 

where F e is the set of equivalent local martingale measures. 

Using the exponential utility for discounted wealth, it is easy to show that 



and substituting the specific expressions of the duals for the buyer and the 
Merton problems, we can write 

P h Q {T)=V--mi Q&f U{Q\P) (28) 



The entropy term in ([28]) can be combined into one entropy term with a 
different prior measure, the minimal entropy martingale measure, which is 
the measure minimizing the relative entropy in Ff. Let 
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g° = arg min Q&f U(Q\P) (29) 
Theorem 2. \Frittell% (200(1) and lDelbaen et al. I \200A) ) Under assumpt 



Q° exists, is unique, is in P/(-P) D P e and its density has the form 

dQ° 



dP 



- 7 x° 



ion 



(30) 



where X® is the optimal terminal wealth associated with the solution of the 
Merton optimization problem (TjJ) and logao = H(Q°\P) < 00 

Under the a ssump tion that ^p- G L 2 (P), we can apply Proposition 2.2 of 



Ilhan et al.l (120081 ) and the bid price is 



Po(T) = inf, 



QePf(Q°) 



E Q {c} + -U(Q\Q C 



1 



(31) 



Next, we define the Radon-Nikodym derivative b etwee n Q° and P for our 



Ilhan et al. 


( 


2008) and 


Bielecki et al. 



(120041 ) . The so-called minimal martingale measure, P°, is defined by the 
following Girsanov transformation. 



dP° 



dP 

The dynamic of (S, Y) under P° are: 



JL 

2a- 



dS t 



rS t dt + aS t dW i 



dY t = (b(Y t )-fa(Y t ))dt + a(Y t )(pdWf oil) +p>dW t Poi2) ), 



where 



w) 



H and Wf o(2) = W t 2 



P° has finite relative entropy and is equivalent to P. Therefore, Q° G P/(P)n 
P fi and is unique. The parameters p and a are independent of Y and by 
Ilhan et al.l (120081 ). the minimal entropy martingale measure coincides with 
the minimal martingale measure. 
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3.5 Comparison with the classical reduced model 

In this section, we compare the indifference price for both buyer and seller 
with the classical reduced price. In fact, the latter is just the limit case of 
our model. Let us recall the classical results of the intensity model. See 



Bielecki and Rutkowskil ( 12002 ) . 

Let Q be a risk neutral probability and the corresponding risk-neutral 
intensity. The price at time t = of the defaultable bond denoted Pq t (T) is 
thus equal to the expectation under Q of the discount pay-off, 



P-(T) = E Q (e- rT l {T>T} ) = e~ rT E Q (e-f° 

Proposition 2. Under the assumption that the market price of risk for the 
non-traded asset Y under Q is equal to zero, 



lim P 6 (T) = lim P°(T) = P~(T) if /, = (33) 

7— >o 7— >o 



Proof. When the risk aversion parameter 7 tends to zero, iBecherer (1200 lh 



proved the indifference price goes to the arbitrage free pricing under the 
minimal entropy martingale measure: 

lim 7 ^ P b (T)=lim 7 ^ P s (T) = E Q ° {e~^ l {r>T} ) 

= E p °(e- rT l {T>T} ) 

If the excess return \i = 0, the minimal martingale measure P° coincides 
with the the investor's subjective measure P and the result follows. 

□ 

3.6 Solving the PDE 

With the change of time t — > T — t, the equation (Tl8l) reduces to: 

-u t + Fi(y)u yy + F 2 (y)u y + F 3 (y)u + F 4 (y)u- e = 
u(0,y) = l 

with : 

Fx{y) = \a{yf 
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F 2 (y) = Ky) - ^a(y) 
a 

FM = - (i - P 2 ) + 
F 4 (y) = (1 - p 2 ) A(y) 

3.6.1 Finite Differencing of PDE 

Our aim is to solve f fT8|) - fl20|) by implicit finite difference schemes. To this 
end, we divide the interval [y m i n , y ma x] into the sub intervals 

y m in = yo <yi— < yjv+i = y ma x 

and we assume for convenience for the mesh-points {yj} N =Q are equidistant 
that is, 

y 3 = + Ay, j = 1, N+l (Ay = v™*-^ ). 

Furthermore, we divide the interval [0, T] into M + 1 equal sub-intervals 

= to < t\... < tu+i = T 

where 

U = t;_! + At, i = 1 , M + 1 (At = jfe) 

Let u(ti + i,yj) = the essence of the finite difference approach lies in 

replacing the derivatives in ffTHj) by divided differences at the mesh-points 
(ti,yj). The following equations are used to approximate the first and sec- 
ond order derivatives: 

8u _ uf 1 - u) 
~dt ~ At 

o i+1 i+1 

att - u j _ 1 

dy ~ 2A~y 

d 2 u _ uf^ - 2^ +1 + u£\ 
W~ (Ay) 2 
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To avoid the difficulty of the nonlinear term, the following first order approx- 
imation is used 

(u) +l y e = (i + e)(u)y e - e{u) +l ){u)y e - 1 . 

Applying the above approximations to ( |34"1) . we get 

a 3 u) + _\ + P 3 u) +1 + v u% = -u) - Uu)Y\ 

with 

AtF x>j AtF 2J 
Ay 2 2Ay ' 

-1 - + AtF 3J - 0AIIYA><;,) " '• 

AtF ld AtF 2J 
Ay 2 2Ay ' 

(i + e)AtF 4d . 

Thus we obtain the linear system Ru l+1 — ef, < i < M where R is 
formed matrix using the coefficients dj, /3j, v,j and d % is a vector with dj = 
—iij — j(u l j)~° . The calculations are similar for w b and w c , the only changes 
are F A (y) = e 7C (1 — p 2 ) \{y) for w b and -F4(?/) = e~ 7C (1 — p 2 ) \{y) for w c . 



a J = 

Pi = 

v i = 

^3 = 



3.6.2 Choice of parameters and Boundaries 

The choice of the appropriate level of risk aversion 7 (> 0) is important 
in the indifference pricing within an exponential utility function. If a high 
value is chosen, the investor will be too risk averse and therefore the results 
of the pricing will converge to the super-replication's prices. On the other 
hand if the value of risk aversion is too small, the decisions of the investor 
will not diff e r app reciably from the decisions of a risk neutral investor. In 



Binici et al.l (120031 ). the authors examine the risk attitudes of farmers in the 
Lower Seyhan Plain of Turkey. While some variation by utility function 
exists in the classification of the sampled farmers into risk averse and risk 
preferring categories, the overwhelming evidence is that the sampled farmers 
are risk averse. The study also found for the exponential utility function, the 
estimated values of 7 G [0.0375; 0.4920] with a mean value of 0.11. We use 
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in our study different values of risk aversion coefficients such as 7 = 0.01, 



7 = 0.2 and 7 = 0.7. We choose as I Sircar and Zariphopouloul (120071 ) p 



0.09, a = 0.15, r = 0.03. We work with negative correlation coefficients 
p G {—0.90, —0.5, —0.20, —0.05} in order to specify that the intensity tends 
to rise when the stock price falls. We define At = Y t with Y t follows CIR 
dynamics defined below. 

Cox-Ingersoll-Ross intensities In this model, the intensity follows the 
stochastic differential equation 



d\ t = a(A - X t )dt + <j>y/\ t dW i 



-3 



where a, A and (f) are positive constants and dWf = pdW} + p'dWj*. When 
we impose the condition 2aA > <p 2 then the intensity A is always positive, 
otherwise we can only guarantee that it is non-negative (with a positive prob- 
ability to terminate to zero). In fact, if the intensity approaches zero then 
the volatility 4>\fXt approaches zero cancelling the effect of the randomness, 
so the intensity rate remains always non-negative. We fix X min = and 



When X t = 1, the agent will follow the sub-optimal policy in investing exclu- 
sively in the default-free bank account (ir t = 0). 

M(t, x, X = 1) = — e _7X and therefore u(t, X max ) = 1. 



H\t,x,X = l) = E(-e-^ +c h {Tt>T} + (-e~^)l {Tt < T} \X t = x,X t = l) 

= -e" 7X + {e~ lx - e-^ x+c) )P{r t >T\X t = 1) 

P{r t >T|A t =l) = l {T>t} £(e-^ A ^|A t = l) 

= l {T>t} A(t,T)e- B ^ 
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with 



B{t,T) 
A(t,T) 

Then for the buyer, 



2( e ^-*) - 1) 
2e+(a + 0(e« T -*)-l) 
2^ e (*K)(^) 

2e + (a + o(^ T -*)-i; 

^a 2 + 20 2 



w b (t, X max ) = (e' c - (e^ c - l)l {T>t} A(t,T)e- B ^)) 1 P * 
Similarly for the seller, 

w'(t, \ max ) = (e^ c - (e^ c - l)l {T>t} A(t,T)e- B ^y- p2 

See Cox et al. fll985l ) for the derivation of A(t, T) and B(t, T). 

When A t = 0, ( IT8l) gives the following equation 



u(T, 0) = 1 



du 



(35) 



The equation at hand is a linear Hamilton- Jacobi equation with the hamil- 
tonian is a function of || and it. The solution of (135]) is 



u(t,0) = e 

Similarly, for the buyer and seller 



w%0)=w s (t,0) = e- {1 - p2) ^ {T - t) 

For the implementation of the model, we use the following estimated values 
for CIR-intensity given in lLongstaff et all ( 120031 ): a = 0.20, A = 0.06, <p = 
0.03. 
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Numerical Results and Comments In this paragraph, we present the 
results of the indifference pricing of the defaultable bond solved by the finite 
difference method. The yield-spread curve which reflects the relation between 
yield spreads and the time to maturity is plotted and analyzed. The yield- 
spread is the difference between the yield on the defaultable bond and the 
yield on the risk free bond (in our situation r). The yield spread is the 
indication of the risk premium required for investing in the risky bond. When 
it is high, this means that the investors require a higher risk premium, so the 
bond is more risky. Conversely when it is small, this means that the bond is 
riskless. 

In the figures H] and [31 we present the spread curves for the buyer and the seller 
respectively for various risk aversion coefficients given the intensity A = 0, 
A = 0.02, A = 0.12 and A = 0.5. When A ^ the short term limit of the 
yield spread for the buyer and seller is nonzero, reflecting the presence of non- 
predictable defaults. Even if the intensity and the risk aversion parameter are 
small, the investor requires at short term a risk premium. We can remark 
for the maturity T = 0, A = 0.02 and 7 = 0.01, the bid and ask spreads 
are equal to 0.045. When the default intensity increases, the bid and ask 
spread curves are in fact high in order to compensate the additional risk; but 
the shape of the curves changes depending if it is a buyer or a seller. An 
experiment on many values of intensity shows the following shapes for the 
buyer's spread curve in the figure [21 

• < A < 0.06, the spread curves are upward sloping reflecting a low 
default risk at short term, whereas the forecast of the credit quality of 
the firm over longer maturities is less certain. This shape is illustrated 
for A = 0.02. 

• 0.06 < A < 0.15, the spread curves are S shaped that is the buyer 
expects that the credit quality of the firm will increase until a cer- 
tain maturity and after that it will be low for longer maturity. This 
particular shape is shown for A = 0.12. 

• 0.15 < A < 1, the spread curves are downward sloping reflecting a 
high default probability of the firm at short term; once the firm has 
survived a certain period of time without a default, it faces a lower 
default probability in the long term. The shape is shown for A = 0.5. 

For the seller's spread curves, we have two kinds of shape in the figure El 
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• < A < 0.25, the curve is upward sloping and is shown for A = 0.02. 

• 0.25 < A < 1, the curves are humped meaning that the seller expects 
to offer a high risk premium until a certain maturity. Once the firm 
doesn't default until this reference maturity, the ask yield-spread de- 
creases. We show the humped curves for A = 0.50 and A = 0.90. 

For A = 0, both buyer and seller's spreads are closed to zero whatever the 
maturity since there is no default probability. The impact of the risk aversion 
7 on the yield-spread is also analyzed. For the buyer, we found a positive re- 
lationship between the spread and the risk aversion coefficient in contrast to 
the negative relationship for the seller. It means that if the buyer is more risk 
averse, he will require a high risk premium for investing in the defaultable 
bond. In contrast, a more risk averse seller will offer a low risk premium for 
selling the defaultable bond. This result is also a consequence of the mono- 
tonicity of the i ndifference p r ice as a function of risk aversion parameter. For 



more details see lllhan et al.l (120081 ). 



We analyze in figures H] and [51 the impact of the correlation coefficient p 
on the spread curves. At the maturity T = 0, the correlation between the 
intensity and the stock price has no impact on the yield-spread for both 
buyer and seller. But for longer maturities, p (in absolute value) increases 
the yield-spreads for the investors. At time T = 50, we can clearly remark 
the impact of p on the bid and ask's spread curves. When \p\ = 1, the 
intensity and the stock price are perfectly linearly dependent and the unique 
source of uncertainty is the stock price which is a traded asset. This leads 
to the complete market. In the figures and [3, we show the convergence of 
the indifference spread to the classical spread for p = and 7 = 0.0001 
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Figure 4: Impact of p on the bid spread-curves 
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Figure 5: Impact of p on the ask spread-curves 




Figure 6: Indifference Spreads and Classical Spreads for buyer, \i 
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Figure 7: Indifference Spreads and Classical Spreads for seller, \i = 0, 7 = 
0.0001 



4 Conclusion 

The classical reduced models for credit risk were modified by embedding 
them in the framework of the utility based valuation. The reason of doing 
this, as mentioned in the introduction, is that the default intensity is a source 
of uncertainty that is not traded and the equivalent martingale measure is 
not unique. This leads to work with the utility indifference valuation for 
defaultable bonds. In our model, besides the uncertainty of the intensity, we 
introduce another correlated diffusion process, which comes from the firm's 
stock price to drive the credit risk. The behavior of the spread curve in 
our model is different from the classical reduced models due to the intro- 
duction of the diffusion process. Two news parameters are introduced in 
the model, namely the risk aversion parameter of the investor and the cor- 
relation coefficient between the intensity and the stock price, which results 
in the nonlinearity of the pricing rule. We derived the Hamilton- Jacobi- 
Bellman equations (HJB) verified by the values functions, and reduced them 
to the nonlinear reaction-diffusion equations. The nonlinear partial differen- 
tial equations were solved by the implicit finite difference scheme. The effects 
of the intensity, the risk aversion parameter and the correlation coefficient on 
the spread curves when analyzed. Although the default intensity increases 
the bid and ask's yield spread, it changes the shape of the spread curves. 
For the buyer, the curves can be upward sloping, downward sloping and S 
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shaped depending of the values for the intensity. For the seller, we got only 
an upward sloping and humped spread curve. The other analysis reveal a 
positive relationship between the bid spread and the risk aversion coefficient 
and a negative relationship between the ask spread and the risk aversion 
coefficient. It reveals that the correlation coefficient increases the bid and 
ask 's spreads only for longer maturities. Another important topic, which is 
not considered in this paper is the asymptotic behavior of the spread curves 
within the stochastic intensity and the pricing of the Credit Default Swap. 
Both of these problems will be analyzed in future work. 
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